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Abstract 
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from the perspective of Bohmian mechanics. A detection probability for detectors sensi- 
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alises an earlier proposal by Leavens and McKinnon. The difference between the two 
notions is illustrated through application to a free wave packet. 
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1 Introduction 

According to quantum theory, the probability of the macroscopic event, which is caused by 
a microscopic system during an act of measurement, is of the type Tr(p to i?). Here E is an 
orthogonal projection within the system's Hilbert space and pt is the system's density operator 
at time to- The time to, at which the approximately instantaneous measurement interaction 
takes place, is determined by the experiment's design. There are, however, important situations 
which do not - even approximately - fit into this framework in any obvious way. Consider for 
instance an unstable nucleus, which is monitored for several days by some initially activated 
detector. You patiently sit next to the detector and register the time, when you hear the click. 
What is the probability that you hear the click during a certain time interval? A cleaner model 
situation involves a freely propagating one particle wave packet, which slowly sweeps over a 
detector activated at time 0. The detector is small compared to the wave packet's size. What is 
the probability P(T) of a click, happening at any time t in the range < t < T? The funcion 
P is called arrival time distribution. Several proposals try to answer this question without 
reproducing the quantum Zeno paradox jjj. Let us describe them briefly. For an extensive 
summary of the subject see §]. 

A first approach attempts to fit P(T) into the scheme Tr(p to E) through quantising the phase 
space function A, which represents the classical time of arrival. The function A maps each phase 
space point from its domain onto that finite time, at which (according to the system's classical 
dynamics) this point enters the detector's location. With E being the spectral projection of 
the quantised A associated with the spectral interval (0, T) it is assumed that P(T) = Tr(poE). 
Working out this general idea reveals that ad hoc regularisation assumptions are needed, in 
order to obtain a self adjoint quantisation of A |||. The need for regularisation is both due to 
the unboundedness of A around p — and to a classically unspecified operator ordering. Due 
to its regularisation ambiguity, this definition of P(T) does not seem convincing. 

A second strategy attempts to derive the arrival time distribution from a unitary quantum 
dynamical model of the continuing observation process and a single final measurement, i.e. 
the "reading out of the observer's notices". To this end, an auxilliary quantum system is 
coupled to the particle during the time interval (0,T). The auxilliary system's position is 
taken as the pointer position of a clock and its evolution is stopped through an interaction 
with the particle's wave function [Q. This approach does not yield the picture of a sudden 
click happening at a certain time t € (0,T), but rather of a smooth influence being exerted 
onto a position distribution. Only the final observation at a controllable time to > T then 
produces the stochastic position value, which is interpreted as an approximate time of arrival. 
Therefore, our macroscopic impression that facts are permanently established in the course of 
time, instead of being created with a final measurement only, - "at the end of the day" , as 
Sheldon Goldstein has phrased itQ - remains unexplained. 

A third way of defining P is obtained by exposing the particle's wave function to an absorb- 
ing detector, whose influence onto the wave function is modelled by a nonhermitean Hamil- 
tonian, e.g. g. Thereby the one particle dynamics becomes nonunitary and the quantity 

— \dt\ for HV'oll = 1 is then - up to an overall normalisation - interpreted as the probabil- 

ity density of clicks at time t. However in general, though not in [||, — may take negative 

values, which in turn implies that the probability P(T) = 1 — H^tII (for the detector to click 
sometimes between time and time T) may decrease upon increasing T. Clearly, a decreasing 
probability is questionable if one imagines say 1000 independent copies of the system side by 
side and the percentage of counters having made their click is monitored as a function of time. 
If this percentage decreases with T, a mechanism seems to be at work, which makes clicks 
unhappened! 

1 private communication 
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Finally, Leavens || , and McKinnon and Leavens |8| have defined an arrival time distribu- 
tion P, which is motivated by the Bohmian flow connected with a solution ip t of Schrodinger's 
equation. They considered the one dimensional case and argued that the (conditional) proba- 
bility density of clicks equals const ■ \ j(t, L)dt\, where j(t, L) is the spatial probability current 
density at the detector's location L at time t. In case of / °° dt \j(t, L)\ —: c < oo, the condi- 



tional arrival time distribution 



T 



P(T) = -- I dt\j(t,L)\ (1) 



is a nondecreasing (nonnegative) function on the interval (0, oo). Yet the integral J °° dt \j(t, L)\ 
need not be finite, as e.g. in the case of an harmonic oscillator dynamics, where the mapping 
t l— * 3 '(*> L) is periodic. In such cases therefore, the definition (0) does not make sense. 

The probability density ~ \j(t,L)\ is derived by Leavens as the "infinitesimal" probability 
that the particle's Bohmian trajectory passes the point L during dt at time t provided the 
Bohmian position at time is distributed by \il>o\ 2 \dx\. If one assumes that the detector clicks 
each time it intersects with the particle's Bohmian trajectory, the density \j(t, L)dt\ j c indeed 
yields the probability density of clicks. This seems to be a reasonable idealisation if the detector 
is active during a short time interval. What happens, however, if the detector is active over a 
longer period of time, such that the same trajectories pass the detector more than once? Does 
a detector really increase its click probability when trajectories intersect, which have done so 
before'? 

The possibility of multiple intersections between Bohmian trajectories and detector posi- 
tions has already been taken into account in [^), JlOfl . These works have ruled out multiple 
intersections for scattering situations. For a summary see also chapter 16 of Diirr's recent 
text book [jllj or As a next step, in the context of "near flield scattering" the exit time 
statistics from a large but finite sphere around the scattering center has been investigated in 
Jt3| . In case of multiple crossings of the sphere's surface by Bohmian trajectories, replacement 
of j in equation ([l]) by a truncated current has been proposed in Jl3| , in order to obtain the 
correct exit time statistics. The truncated current only counts the first exit of trajectories as 
detection events. This is reasonable if the initial wave packet is well localised within the sphere. 
Yet if a considerable part of the wave packet has left the sphere by the time the detector is 
activated, those trajectories, which have already entered the detector and stay there, carry 
a nonncgligible portion of probability. Accordingly they should contribute to the detector's 
click probability. However these trajectories do not contribute to the surface integral of the 
truncated current. Thus under such circumstances a more general prescription is needed to 
count also those trajectories which are confined to the detection volume during the full period 
of detector activity. 

In this work we propose and explore a very natural definition of detection probability within 
Bohmian mechanics, which on the one hand meets the above needs for generalisation but also 
implies the idea of using the truncated current for the exit time problem described in |]l3| . 
The physical argument behind it is quite simple: a realised trajectory induces a detection 
event at the earliest instance only, when this trajectory falls into the detector's volume because 
thereafter the detector remains discharged. According to Bohmian mechanics each individual 
trajectory is realised with an "infinitesimal" probability to be computed from the wave function. 
Adding up these infinitesimal probabilities for all the trajectories intersecting the detector's 
volume during its period of activity then yields this detector's click probability. Assuming 
this, we obtain an expression for the arrival time probability density, which in general depends 
on the spacetime region to which the detector is sensitive. In cases, where each Bohmian 
trajectory crosses a (point like) detector at most once, equation (Q) remains valid. However, 
our definition yields different probabilities otherwise. 

Why an experimental decision between the various conflicting proposals for P(T) is not 
yet feasible has been indicated in section 10 of the latest review of the subject ||. The basic 
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reason seems to lie in the difficulties in preparing a specific wave packet which has to be large 
compared to the detector's size and which in addition has to pass the detector sufficiently 
slowly. Clearly such experiments are not precluded on principle. 

After a brief summary of Bohmian mechanics in section 2, we develop our definition of the 
arrival time distribution P(T) in sections 3, 4, and 5 within a Galilean spacetime frame work. 
Using spacetime proves very suggestive since Bohmian trajectories become one dimensional 
submanifolds (worldlines) instead of mappings. In section 3 we define the Galilean (one particle) 
spacetime from its structural atlas. Section 4 contains an outline of conserved flows on Galilean 
spacetime. Here again we choose the coordinate independent formalism of exterior calculus. 
This has the following advantage. When computing the flux through a (possibly moving) 
hypersurface it is the current 3-form which is integrated over a 3-manifold. Neither a metric 
spacetime structure nor a normal vector field, both breaking Galilean invariance, need to be 
introduced. The main result of this section is the formula of definition 4 which gives the 
amount of conserved "mass" passing through a spacetime region X. In section 5 we apply 
this formula to the flow of the quantum mechanical position probability. Here our definition 
of the detection probability is obtained from the quantum mechanical probability measure on 
the set of Bohmian worldlines (orbit space), which follows from the wave function ip under 
consideration. The^ probability that a detector clicks, is assumed to equal the probability 
measure of the set X of all those orbits, that have a nonempty intersection with the spacetime 
region X , to which the detector is sensitive. This measure in turn equals the usual quantum 
mechanical probability measure to detect a particle with wave function ip at t = within the 
set of all those locations which are taken by the orbits of X at t = 0. Our definition works 
for very general, extended spacetime regions and it works for the free Schrodinger dynamics 
as well as for ones with nonzero potential. The coordinate independent treatment guarantees 
Galilean invariance of the detection probability for zero potential. In section 6 we illustrate 
our notion of P(T) through the example of a free standing Gaussian wave packet. 

Our definition of P(T) in terms of the Bohmian flow could be improved by taking into 
account the detector's influence onto the particle's Bohmian trajectories. Since the latter 
become projections of the higher dimensional orbits of the detector plus particle system, this 
effect can be considerable even for detectors without any back reaction onto the particle wave 
function |fl4| . The general idea of our approach however, remains the same. Also if the detector 
(or even the observer) is modeled as part of the quantum system, an assumption has to be 
made about when each individual orbit generates the click (in the observer's mind). This rule 
then mathematically represents the discrete event, which is missing from standard quantum 
theory. 

2 Summary of Bohmian Mechanics 

The density operator p, representing the state of a quantum system with (separable) Hilbert 
space H, defines a probability measure W Pt A on the spectrum of any self adjoint Operator 
A of TL. It is given by W p _a(P) — Tt(pEa(I)) where Ea{I) denotes the spectral projection 
of A associated with the Borel set I C spec(A) C M. Standard quantum theory assumes 
that, if a measurement of the observable A is performed on the state p, then W p ,a{I) equals 
the probability of the event "the measured spectral value a of A belongs to J". Now, for 
dim(K) > 2, there does not exist a density operator p such that W p ^a is dispersion free, i.e. 
a point measure, for all A. Gleason |L5| has investigated the question whether there exist 
more general ways of defining a probability measure on spec (A) for all A. To this end he 
considered the mappings a from the set n of all orthogonal projections of 7i into the real 
numbers such that <r(P) > for all P e n and cr(J2i Pi) = a (Pi) I0r am/ countable sum of 
Pi 6 n, with PiPj — Si.jPi. In addition he assumed a(id) — 1. From this he derived in case of 
dim(7i) > 3 that for any such mapping a : II — > [0, 1] there exists a density operator p such that 
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a{P) = Tr(pP) for all Pen. Thus the idea of generalising the formula W pA {I) = Tv(pE A {I)) 
to Wct.a(I) — <t (Pa(^)) in order to possibly obtain "deterministic states", i.e. point measures 
W„ t A f° r au A, and a representation of density operators as mixtures of these, under the 
adopted assumptions fails. 

The standard quantum physical interpretation of this body of mathematical facts leads to 
the following conclusion. It is inconsistent to suppose that the state of an individual quantum 
system is a deterministic state, i.e. determines values for all observables, and it is inconsistent 
to suppose that a density operator p only describes a mixture of such fictitious deterministic 
states. (It is generally held inconsistent to suppose that an idividual particle has a specific 
position and a specific momentum and so on.) 

From this conclusion then the notorious quantum measurement problem follows: How can 
standard quantum theory represent within its formalism the mere fact that individual closed 
systems do have properties? (This surely is the case for systems comprising an observer and 
not beeing in need of any sort of external observation inducing a state reduction, the quantised 
deus ex machina.) Which fact concerning a closed system is it, whose probability of being the 
case is given by W p ,a{I) • 

Bohmian mechanics resolves these problems for systems with a Schrodinger (or Dirac) 
equation: a picture of individual systems with defined properties emerges. A concise review of 
Bohmian mechanics is given in reference [ Jul . An informal but clear summary is to be found 
in Goldstein's contribution to the Stanford Encyclopedia of Philosophy jL7| . Let us summarise 
the basic ideas. 

Bohmian mechanics introduces deterministic states which violate Gleason's assumptions 
and accordingly circumvent his theorem. It is assumed that an individual system has a state 
(ip, q) given by a wave function t/j in the system's Hilbert space and a point q in its config- 
uration space, q is supposed to represent the actual positions of the system's constituents. 
Other observable properties of the system have to be derived from the Bohmian state through 
a dynamical analysis of the concrete experiment designed to measure them. In this way all 
other properties like spin or momentum are expressed through the state's well defined posi- 
tion properties. It turns out, that the spectral value, which a general observable assumes in a 
Bohmian state, depends on the specific way of how this observable is measured, i.e. contex- 
tuality is found to be realised [fl8|| . Accordingly Gleason's assumptions on the mapping a are 
violated because a needs a much more complex domain than simply the set of all orthogonal 
projections. 

In order to work out the dynamical program of reducing all state properties to position 
properties, the time evolution of Bohmian states is needed. It is assumed to be given by a 
Schrodinger equation for the wave function and by a time dependent tangent vector field v on 
the configuration space, v is defined in terms of the solution ip t of the adopted Schrodinger 
equation with initial condition ipQ — ip. The integral curve j q of v with initial condition 
7g(0) = q gives the system's configuration at time t by j q (t). 

Finally, Bohmian mechanics establishes contact with empirical data. This happens accord- 
ing to the rule of quantum equilibrium. It states that for an ensemble of systems, each with 
wave function ip, the individual system's position q belongs to a configuration space domain 
A with the usual probability J A \ip\ 2 d n q. A controlled preparation of q contradicting this rule 
is assumed to be impossible by present day technology. (All this can be justified to a certain 
extent within the Bohmian picture [jl6| . ) And finally as a last ingredient it is supposed that it 
is the center of mass position of pointers and the like that we observe. 

The rules of Bohmian mechanics are such that the probabilistic statements of standard 
quantum theory are reproduced. So there seems no room left to argue about the empirical 
superiority of either standard quantum mechanics or its Bohmian extension at the ensemble 
level. Bohmian mechanics might, however, give a clue for the correct treatment of ensemble 
problems, where the standard interpretation remains unclear and offers conflicting strategies. 
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As described above, standard quantum theory offers various different conceptions for the arrival 
time distribution P(T). Therefore we hold the arrival time problem to be one such opportunity 
for Bohmian mechanics to possibly show that it also has its value in dealing with ensemble 
problems on top of its merit of providing a language for speaking about individual systems. 
We add another conception for P(T) which is motivated by the Bohmian extension of quantum 
mechanics. It does not conform to the standard scheme of identifying P(T) with some quantity 
of the type Tv(pE) with E being independent from p. Our P(T) needs the concept of Bohmian 
trajectories for its very formulation. One should note, however, that Bohmian trajectories are 
implicit in the wave function, whether one intends to make use of them or not. 

3 Galilean Spacetime 

We model spacetime as a Galilean manifold. Various equivalent definitions of a Galilean man- 
ifold can be given. Here we use the method of a structural atlas. The basic object is the group 
Q of (orthochronous) Galilei transformations. 



(2) 



r := eGi„+i(R) \veR n ,ReO n y 

g := { 9 : ^ R ,i+1 , = 7 ^ + a | 7 e T, a e R" +1 } . 

The elements of R™ +1 and R" are treated as column vectors throughout the text. 

Definition 1 A Galilean manifold (M,Ag) consists of a differ entiable manifold M. and a 
subset Ag of the atlas A of M., where Ag contains global charts only and the set of transition 
functions {$2 | 3>i,3>2 € Ag} equals Q. The charts $ G Ag are called Galilean charts. 

A Galilean manifold carries the canonical time-l-form 9 := with $ = (<£>°, $ 1 , ..<&")* € 
Ag. Observe that 6 is independent from the choice of Tangent vectors v € T{M) with 
9{v) = 1 are called velocity vectors, and tangent vectors with 6(v) — are called spacelike 
vectors. The subbundle TZ(Ai) := ker(#) of spacelike vectors is completely integrable. Its 
integral manifolds are given by S$ ;t := {p e M | $°(f>) = t) where $ E Ag and t £ R. These 
integral manifolds are called instantaneous spaces. 

The vector bundle 1Z(M) carries a canonical positive definite fibre metric 



n 

(;•) :=V d<$> k ®d$ k , 



fe=i 



where the restriction of d<fr k to 1Z{M) is again denoted as d$ fe . Note that (•, •) is well defined 
as a fibre metric of 1Z(M), but is not so as a fibre metric of T(M). Finally, the Galilean 
manifold carries two orientations represented by the two volume (n + l)-forms 

Q := {±d$° A d® 1 A ... A d^ n } , $ G Ag ■ 

Thus the density \d$° A dQ 1 A ... A is unique. Various further structures are canonically 
defined on (M,Ag) as e.g. a linear connection of the tangent bundle. We shall not use them 
here. 

4 Conserved Flows 

Let j be a diffcrcntiablc tangent vector field on a Galilean manifold (Ai,Ag). By choosing a 
volume form to G £1 the differentiable n-form J on M is obtained through 

J :=Jjuj : (ii,.. t n ) h-> u(j,ti, ..t n ). 
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The associated density |J| does not depend on the chosen u>. The divergence of j is the unique 
function (see e.g. page 281 of div(j) satisfying 

LjLu = div(j)cj. (3) 

Here LjU denotes the Lie derivative of uj with respect to j. The divergence of j does not 
depend on the choice w e fl. Observe that this definition of the divergence of a vector field 
does not make use of a (pseudo-) Riemannian metric. It is built on a given density As 
it is the case with a Galilean manifold, this density \ui\ need not be induced by a (pseudo-) 
Riemannian metric. If |w| is the metric density of a (pseudo-) Riemannian manifold, the above 
definition for div(j) coincides with the usual one. 

There holds Ljtu — jjdeo + d(jjco) — dJ , and therefore 

div(j) = ^ dJ = 0. 

Furthermore, if div(j) = 0, then LjJ = j-idJ + d(jjJ) = j-idJ = 0. In terms of a Galilean 
chart $ the divergence of the vector field j = J2k=o 3% ' ®t rea ds 

n 

div(i) = E 9 *b'£]- 

fc=0 

Here 9* denotes the tangent vector field associated with the chart <& and are the 
coefficient functions of j with respect to the coordinate frame 9* :— (9*,. ..9*). For u> = 
d&° A dQ 1 A ... A d<& n the n-form J is given by 

n 

J = ^(-1) & • j% ■ A ... A d^- 1 A d$ fc+1 ... A d$ n . 

Let j be a C 1 -vector field on M. such that p (j) ^ for all p e M. Then the velocity 
vector field of j is defined on M by j := ~ghxj- The maximal integral curve of j through p £ M 
is the (unique) function 7 : I — > M. with 7(0) — p and 

T(A) = j 7 (A) for all 

Here / is an open real interval, which cannot be extended. The image 7(1) C M. is called 
(integral) orbit of j through p. 

Assume the vector field j on M to be complete, i.e. each maximal integral curve of j 
has the domain M. Then a unique one parameter group of mappings F s : M. — * M with 
F s (p) — 7(s) exists, where 7 is the maximal integral curve of j with 7(0) = p. There holds 
F s oF t = F s+t and Ff 1 = F_ s for all s.tel. The mapping F:lxM^ M, (s,p) h-> F s (p) 
is called the flow of j. Since 6(J) — 1, for the maximal integral curve through any p G M. there 
holds ($° o 7) (s) = &°(p) + s for any s € M and for any $ S Thus no orbit begins or 
ends at finite time. In particular F s carries instantaneous spaces into instantaneous spaces, i.e. 

If now div(j) = 0, we have dJ = 0. From this and because of j_iJ = j_i(j_iuj) — there 
follows LjJ — j-idJ + d(gjjjjjJ) = and therefore both J and | J| are invariant under the pull 

back with the flow of j, i.e. F*J — J and also F* \ J\ — \ J\ for all s £ R. From the pull back 
formula for integrals of differential forms then the following lemma follows. 

Lemma 2 (Integral conservation law) Let j be a C 1 -vector field on M. such that div(j) = 
0, 9 p {j) ^ for all p £ M and such that j is complete. F denote the flow of j . Then for any 
Borel set of an instantaneous space X C and for any i£ i there holds 

I \A= I \J\- (4) 

JF t (X) JX 
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Remark 3 Depending on the physical context an integral of the type J x \J\ is interpreted as 
the mass or probability "contained" in the instantaneous region X. The lemma thus establishes 
the picture of a flow which transports mass or probability without change along the flow lines. 
The same amount of mass which is contained in an instantaneous region X is contained in 
F t (X) for any t G E. 

Consider now more general sets X C M. which need not be contained in an instantaneous 
subspace. Let us try to formulate a precise notion of the amount of mass passing through X. 
A clear and unambiguous way of doing this is by determining the set X of all orbits passing 
through X and by computing the amount of mass carried by these orbits. This can be done 
by intersecting these orbits with any instantaneous space E$ i4 and by integrating \J\ over this 
intersection. Thus we have motivated the following definition, which is illustrated by figure 1. 

Definition 4 Let j be a C 1 -vector field on M. such that div(j) = 0, p (j) ^ for all p £ M. 
and such that j is complete. F denote the flow of j. Let pr be the projection pr : R x A4 — > 
M,(t,p) i — ► p and let E^^ t be the restriction of the flow F foRx£$ jt . Thenir^t '■— P ro is 
the fibre projection of M. onto 2$,* along the orbits of j . If for a subset X of M. its projection 
7r$ jt (X) C is a Borel set, then we define the transition P [X] of j through X as 



Remark 5 Note that the transition P[X] does not depend on the chosen hypersurface 
This follows immediately from ir$ iS+t — Ft ° and from equation fy) because of 



Remark 6 Let X\,X^ G M be disjoint. Then the sets n$j(Xi) and 7r$ j t(X2) need not be 
disjoint. As a consequence P [X± U X2] 7^ P [Xi] +P [X2] in general. Thus P is not a measure. 
Yet Xi C X 2 implies P [Xj] < P [X 2 ]. 






Figure 1: transition of j through X 
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5 Detection probability from Bohmian flow 

In order to define a (free) Schrodinger equation on a Galilean manifold (A4,Ag), one has to 
choose a tangent frame d , which is associated with a Galilean chart <&. Any two such Galilean 
charts $i and $2 are connected by $2 = 9 $1 = 7 • $1 + a with 7 e T and a € R n+1 . The 
frames then obey 9* 2 = 9* 1 • 7. In terms of this matrix notation the duality between a 
frame and its co-frame d<E> := (d$°, ..d^™)' is expressed by the equation d<I>(<9*) = I n +i, with 
I n +i G GZ n+ i(R) being the unit matrix. There holds d<E> 2 = 7 _1 ■ Note that d$ 2 = 
and d* 2 = 9* 1 for 7 = 7„+i, such that a chosen frame determines the chart <f> g .4g up to an 
element a g R n+1 . 

For every Galilean frame <9* we define the differential operator Dg* operating on C 2 - 
functions tp : M. —* C through 



fc2 « 

fc— 1 



The operators D d <s> depend on the frame <9* because of the term ihd® . If $1 and $ 2 ar e two 

Galilean frames with 9* 2 = 9* 1 • 7 and 7 = ^) then there holds <9* 2 = 9* 1 + X^fe=i t,fe ^fc 1 - 
The following proposition however, which can be checked easily, shows that the solution spaces 
ker (.Dg*) can be mapped bijectively onto each other. 

Proposition 7 Let 9 $1 ,9* 2 fee Galilean frames with 9* 2 = 9* 1 • 7, and 7 = (* £et the 

function <p ■ A4 — > R &e given 6y 



ft I 2 



$?-^V$M+c, cel. 



fe=i 



T/ien ker(D a 4. 1 ) is mapped bijectively onto ker(Z? 9 * 2 ) through tp \— > exp(i<fi)ip. 

Let ^ solve the free Schrodinger equation D d i>ip = 0. Then the current (vector field) 
j(V')^*) is defined by 

mJ fc=i 

Due to D d ®tp = there holds div (^j(ip,d®)j = 0. The current's frame independence follows 
through a straight forward computation. 

Proposition 8 Let 9 $1 ,9* 2 be Galilean frames and let ip g ker(D a # 1 ). Then j(tp,d' s ' 1 ) = 
j(exp(i(p)tp,d^ 2 ). 



Remark 9 For Tp g ker(.D a *) we t/ius abbreviate j := j{ip,d®). 

For ^ g ker(.D a *) the unitarity of the Schrodinger evolution implies that the integral 
is* j (V^VO ' A .. A d<&™| is independent of t. If this integral is finite, it may be assumed 
to be equal to 1 without loss of generality. In this case each of the hypersurfaces S$ it carries 
the probability measure defined for the Borel sets X C S$ it 



M t (X) := [ (ip*ip) ■ \d& A .. A d<&™| = f \J\, 
Jx Jx 
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where J = jjco with li chosen from {±d$° A .. A d$™}. The form J is closed because of 
div(i)=0. 

In case of i/j*^ > the vector held j is dehncd on all of A4. If j is complete, its global how 
F provides a hbration of M. by its orbits. The mappings F t evolve instantaneous regions from 
E$. s into instantaneous regions from E$ iS+t of the same probability content. Thus the orbit 
space carries the unique probability measure, given by 

fi(Y) := M t {{x € E*, t | 3o G Y with x G o}) 

for any t G R. Thus for the transition of j through a set X c .M there holds P [X] = G 
[0, 1]. Here X denotes the set of F-orbits intersecting X. 

Bohmian mechanics proposes to take serious the flow lines, i.e. the orbits of j, as the 
possible worldlincs of a quantum point particle with the wave function tp. Which orbit is 
realised in each individual case of an ensemble, is considered as being beyond experimental 
control, and is assumed to be subject to the probability measure represented by Mo. In this 
way Bohmian mechanis provides a picture of a world with facts, evolving continuously in time, 
while simultaneously the quantum mechanical expectation values of fixed time measurements 
remain unaltered. A generalisation of Bohmian mechanics to wave functions, that do not yield 
a globally defined complete velocity vector field, has been established in [Q . 

Within the Bohmian extension of quantum mechanics, the following notion of detection 
probability seems plausible. The probability that the Bohmian orbit of a (free) particle with 
wave function ip G ker(£)g«) passes a given spacetime region X C M, equals the transition 
P [X] of the current vector held j(V , ,<?*) through X. Observe that P [X] does not depend on 
the choice of $ G Ag and that indeed < P [X] < 1 holds. We now suggest that an (idealised) 
detector, which is sensitive to the spacetime region X, registers the particle if and only if the 
particle's Bohmian trajectory passes X. Therefore we assume the detection probability within 
the spacetime region X to equal P [X] . 

Let us consider a more specific situation. Let the set X C M. be the union of time translates 
of a Borel subset D of the instantaneous space S^ o, i-e. 

$(X) := {(t,xf \T 1 <t<T 2 and (0,2:)* G $(£>)} 

for given T\ < T 2 G R. The set X contains the spacetime points covered by an inertial, rigid 
detector, which is activated at time T\ and which is turned off at time T 2 . The number P [X] 
is the probability that this detector clicks. 
The mapping 

5 : {(T lt T 2 ) G K x R | Ti < T 2 } —> [0, 1] , {T X ,T 2 ) h- P [X] 

is continuous. Furthermore the function T 2 i— > S(Ti,T 2 ) is nondecreasing and the function 
T% i ► 5{Ti,T 2 ) is nonincreasing. Thus turning off later with T\ being kept fixed does not 
diminish and activating later with T 2 being kept fixed does not increase the detection probablity. 

In the next section we shall make use of the h — 1 and m — 1 simplihcation of Schrodinger's 
equation. This is obtained by introducing the affine (non Galilean) chart % — (x°iX 1 ^ ■■X n ) — 
(±$°, ...^$ n ). Therefore we have 

d X ° =idd>°, dx 1 = ^d$ 1 ,...dx n = ^d$ n , 
q$ _ lax a* _ Vrnp,X a* _ \Angy 

Then -0 G ker(Z? a *) is equivalent to 

n 
k—l 
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The current vector field j, given by equation (||), and the volume form lu := d<£>° A .. A d$ n 
have the following coordinate expressions in terms of X - 

uj = ^-d X ° A..Adx n - 

Thus in terms of the rescaled wave function := ^^7=) ip the current form J = jjui finally 
reads as follows 

J = ^Vdx 1 A ... A d X n - -J- [** (d?V) - cc] d X ° A d X 2 A .. A d X n + ... 



6 P{T) for a Gaussian wave packet 
6.1 The flow map 

We assume n = 1 in what follows and we use the more suggestive notation: x° =: r an d 
X 1 = : £. Accordingly we abbreviate: <9g = <9 T and — d^. Let 5 G R>o- Then the complex 
valued function ip on A4 



TTl .it 1 

* with * := 

h 



• exp 



il 1 



solves the Schrodinger equation, i.e. Dq^i/j = 0. It is a Gaussian wave packet centered at £ = 
at all times. The complex square root has its cut along the negative real axis. The current 
vector field j := j(ip,d ) is given by 



1 



^ + ^2 9 « 



with 



,7A ex P("!^ 



Here the positive realvalued function A, defined on M, is given by 



A:= ' V1+ (J 



For later use we introduce the rescaled current s := -^=j = s°d T + s 1 d^. The velocity vector 
field associated with j 



j = d T + 



is of C°°-type on M.. 

The integral curves 7 P of the velocity vector field j through a point p G M. are obtained 
in terms of the functions x° :— t o j p and x 1 := £ o -y p . They solve the system of first order 
differential equations 
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with the initial condition p° := x°(0) = r(p) and p 1 := x 1 (Q) — The first differential 

equation has the unique, maximal solution x (A) = X+p° for any A € R. Inserting this solution 
into the second equation yields the non autonomous first order differential equation 



^(A) = 



(p° + X)-x 1 (X) 



J 4 1 



P°+A 
S 2 



Its unique, maximal solution is obtained by separation of variables. It is given by 



x 1 (X)=p 1 



\ 



for any A € R. Thus the vector field j is complete and the Sow F : R x A4 — ► M. defines a one 
parameter group of global diffeomorphisms \F\ \ X £ R} of M. . The coordinate expression of 
F\ is as follows. 



$ o F A o $ 



/ 



2 ,(P°,P 1 )' 



V 



\ 



p«+\\ 2 \ 



J 



The (maximal) integral orbit of j through p £ £$,o is the set of points T p C M on which 
holds £ = p 1 . It is the well known hyperbolic worldlinc of the Bohmian particle with wave 
function ip and passing through p. See e.g. sect. 4. 7 of ref. pl| . Some orbits are shown by 
figure 2 in terms of the dimensionless coordinates t :— t/5 2 and x := £/6. 




Figure 2: Bohmian orbits of a Gaussian wave packet 



The 1-form J = j_iu> = s°d^ — s dr obeys 



J = 



5 2 A- 



dr 



Due to Poincare's lemma, J is exact, i.e. there exist functions H : Ai 
any two functions Hi and H2 with dil\ = di?2 = J the difference i?i - 



-> R with J = rfiJ. For 
i?2 is constant on A4. 
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Due to dH — (d T H) dr + (d^H) d£, for the function H there holds 

d^H = J(d 6 ) = s° = and 

d T H = J{d T ) = s 1 = -^2*** ■ 
A solution to these equations is given by 

2 \A 

where erf : R — * (—1, 1) denotes Gauss's error function 

2 f x 

erf(x) := — = / cxp(— z 2 )dz . 
V 71 " Jo 

Obviously, H is constant on the orbits of j. This is due to dH(j) = J(j) = u>(j,j) = 0. 
6.2 Detector activated at time 

Now we shall discuss the detection probability of a pointlike detector, wich is exposed to the 
wave function ip. The detector is assumed to be located at £ = L > and is activated at r = 0. 
Thus the detector measures the transition of the current j through the spacetime regions 

D T -={peM \ f (p) = L and < r(p) < T} with T > 0. 

The boundary of D T equals {A,B} with {t,£){A) = (0, L) and (r,£)(B) = (T,L) (see figure 
3). 




Figure 3: Detector at rest from A to B 
The set of points po S whose integral orbits r po intersect Dt is, due to A(po) = 

7T (D T ) = lpo e S$ o I there exists a p e D T with |^rr = ^r^l ■ 
I A (p) <5 J 

Thus we obtain 

LS „, . L<S 



7T (D T ) = \p G E$, | -— - < C(po) < 



A(B) " ^ ' " A(A) 
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Due to A(A) = S and A(B) = S\/ 1 + ( j) , this yields 



TT (D T ) = < po G S$ )0 



< e(po) < £ 



1 



The boundary of the line segment tt (Dt) equals {A, C} with 

L 



(r,O(G) = (0, 



?)• 



V 1 + (t) S 

The detection probability P [£>t] then follows by integrating |J| over tt(Dt)- 



- / 






Jtt(Dt) 




1 








erf 


(f 


2 








<5\/l + (?) 2 



=:«l(0,T) 



The function (5^(0,-) is monotonically increasing, has the value at T = and tends to 
|erf(j) S (0, |) for T — > oo. The detection probability stays below 1/2 because no left 
moving orbit intersects with the detection region Dt- The limit of a far away detector yields 
lim L ^oo limbec 5 L {0,T) = 1/2. 

Figure 4 shows P [Dt] as a function of the dimensionless time t := w for i = 1005, i.e. the 
function 

/ : R> -> [0, 1] , i ^> 1 ^erf(lOO) - erf (^=p 




Figure 4: Detction probability P [Dt] 

The equality between P [Dt] and P(T) as given by Leavens @, §, @, we denote it P L (T), 
can be derived as follows. The line segment Dt has the boundary points A and B. The points 
B and C belong to the same orbit Tb of j. The part of Tb lying inbetween B and C is denoted 
by Tb.c Thus the union of the three segments Dt, Tb,c and tt (Dt) is a closed line K C A1. 
The orientation of K and its boundary dK is determined by the chosen u> |p2| . Application of 
Stoke's theorem to the spacetime region K interior to this closed line gives 




/ J+ [ J+ [ J. 

J-k(D t ) Jd t Jt bc 
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Since J r J — 0, because of j j J — 0, and s°, s 1 > on we obtain from this 



P[Dt}= / |J|= / |J|= / l^drl =: P £ (T). 

Jtt(Dt) JD t J D t 

Due to H(B) = H(C), one explicitly verifies 

P £ (T)= / |J|=/ |dff|=/ |(a r fl)dr|=H(A)-H(B)=P[£) r ]. 
Jd t jd t Jd t 

Thus in the present case the detection probability P [Dt] is obtained by integrating the density 
sMt along the detector worldline Dt- Obviously, the equation 



P[D T \ = / [s'dr] (6) 

J D T 

is due to the absence of multiple intersections between Dt and the individual Bohmian orbits. 
We shall construct an explicit counterexample to equation (^) in the next subsection. 

From the function (5^(0, •), the conditional probability density of arrival times at a detector, 
which is activated at t = 0, can be obtained as follows. The conditioning is with respect to those 
events, where the particle is detected at all by this detector. Define the normalised conditional 
distribution function W(T) := ^^j^Ws = "^h ~{A) B) • The differential dW yields the 
conditional probability density w \dT\ := \dW\ of detection times. Thus w(T) = dw ^ 



w(T) 



dT 

(d T H)(B) j\B) 



H{A) H(A) 
' ( -a T ed(^)){U) 



erf(f)V 

2 LT 



/5Ferf(f; 



3 eX P I , 9\ I ' 

1 + \ Hfr)')) 



The density w of the dimensionless time t := T/8 2 is defined through w{t)dt = w(T)dT 
and thus with A := L/S we obtain 

2A t f\ 2 
w(t) = —= — — — — • 3- • exp 



/7rerf(A) (l + t 2 )^ V 1 + * 2 

Figure 5 shows the graph of w for A = 100. 

Since lim^oo t 2 w(t) > 0, the improper integral limA^oo f tw(t)dt does not exist. Thus 
an average (conditional) detection time does not exist as well. 



6.3 Detector activated before time 

In order to be sensitive to the contractive phase of the wave function, we now assume that the 
detector is turned on at some time T4 < 0. It thus measures the transition through the sets 
of spacetime points 

D T := {p e M \ £{p) = L and T A < r(p) < T} with T > T A . 

The bounary dD T equals {A, B}, where (r, g)(A) = (T A , L) with T A < 0, L > and (r, £)(B) = 
(T, L) (see figure 6). We shall see the difference between P [Dt] and P(T) according to Leavens 
0, we again denote it as Pl(T), clearly. 
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Figure 5: Conditional probability density w of arrival times 




Figure 6: Detector at rest from A to B 



The transition P[Dt] then follows by inspection of 7t$ j0 (-Dt)- With the auxilliary point 
C := -DtHE^o we obtain in terms of the dimensionless coordinates t = T/5 2 ,t A = T A /8 2 , A = 
L/6 



P\Dn 



H{B) - H(A) for T A < T < 
H(C) - H(A) for < T < -T A 
H(C) - H(B) for -T A < T 



rf 



-erf 



erf W- erf (7*if)) 

| (erf (A)- erf (-^)) 



for t A < t < 

for < t < -tA 
for — t A < t 



Figure 7 shows P [Dt] (solid line) as a function of t for A = 100 and t A = — Vo • 100. For t > 
our expression P [Dt] for the detection probability P(T) differs considerably from the integral 
of | J | over Dt, proposed by Leavens to represent P{T). This latter integral yields 



Pl(T) 



H(B) - H(A) for t < 

2H{C) - H{A) - H(B) for t > 
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for t A < t < 
for t > 



Its dependence of t is shown for A = 100 and Ia = — V% ■ 100 as a dashed line in figure7. 
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Figure 7: Distribution functions P [Dx] and Pl{T) 

P [Dt] is constant for < T < —Ta, while Pl has a point of stationarity only for T = 0. 
For < T < —Ta orbits cross the detector's worldline, which have done so before. Only past 
the point A 1 with £(A') — L and t(A') = —Ta the probability P [Dt] increases again, because 
orbits are passing, which have not done so before. 

Figure 8 finally shows the conditional probabilities 

P{T) 
limr^ P(T) ' 

associated with Leavens' proposal P(T) = Pl(T) (dashed) and P(T) = P [Dt] (solid) respec- 
tively. 
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Figure 8: Conditional distribution functions of Pl{T) and P [Dt 
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